PL-TR-92-2262 

Environmental  Research  Papers,  No.  1110 


AD-A263  572 


WIND  MEASUREMENT  ACCURACY  OF  A  LORAN 
RADIOSONDE 


James  F.  Morrissey 


26  October  1992 


DTI 

i»JBCT 


>-> 


fllicrr 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


PHILLIPS  LABORATORY 

Directorate  of  Geophysics 

AIR  FORCE  MATERIEL  COMMAND 

HANSCOM  AIR  FORCE  BASE,  MA  01731-5000 


93-08714 

■llllllll 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


W.  SNOW 

fief,  Satellite  Meteorology  Branch 
Atmospheric  Sciences  Division 


Director,  Atmospheric  Sciences  Divisi 


This  report  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center.  All  others  should  apply  to  the  National  Technical  Information 
Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization,  please  notify  PL/TSI, 
Hanscom  AFB,  MA  01731-5000.  This  will  assist  us  in  maintaining  a  current  mailing 
list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0168 

Pubhc  reporting  for  thl«  coV*dton  of  information  la  asti  mated  to  average  1  hour  per  reeponee.  Including  the  time  for  reviewing  inetruetfone,  •#  arching  exit  ting  data  eourcee, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  co  Bee  bon  of  information.  Send  oommente  regarding  this  burden  estimate  or  any  other  aspect  of  this 
colection  of  information,  including  suggestions  for  reducing  this  burdsn,  to  Washington  Hsadquartsrs  Servfcss.  Directorate  for  Information  Operations  and  Reporta.  1215  Jefferson 

Davis  Highway.  Suita  1204,  Arlington,  VA  22202*4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-01 85),  Washington.  DC  20503. 

t  AGENCY  USE  ONLY  (L«n  bl«r*)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

26  October  1992  Scientific  Interim 

4.  TITLE  AND  SUBTITLE 

Wind  Measurement  Accuracy  of  a  Loran  Radiosonde 

5.  FUNDING  NUMBERS 

PE:  62101F 

66701905 

5.  AUTHOR(S) 

James  F.  Morrissey 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSEES) 

Phillips  Laboratory  (GPAS) 

Hanscom  AFB,  MA  01731-5000 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

PL-TR-92-2262 

ERP.  No.  1110 

9  SPONSORING'MONITORINQ  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

1 1 .  SUPPLEMENTARY  NOTES 

12a.  DISTRIBim ON/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUTION  COOE 

13.  ABSTRACT  (Maximum  200  word.) 

A  comparison  of  winds  obtained  using  a  Loran  radiosonde  system  developed  at  Phillips  Laboratory 
and  winds  derived  using  radar  tracking  of  targets  attached  to  the  radiosonde  balloon  train  are 
presented.  The  tests  Involve  10  ascents  with  dual  Loran  ground  stations  on  most  flights  and  dual  radar 
tracking  on  two  flights.  The  tests  wer  performed  at  the  Goddard  Space  Flight  Center /Wallops  Flight 
Facility  In  June  1990.  The  radiosonde  Is  sufficiently  accurate  (0.32  m/s)  to  evaluate  a  wind  profiler  that 
has  an  error  budget  of  1 .0  m/s.  The  wind  profiler  has  been  proposed  to  satisfy  Day  of  Launch  missile 
range  wind  measurement  requirements. 

14  SUBJECT  TERMS 

Wind  measurement.  Radiosondes.  Atmospheric  wind 

15.  NUMBER  OF  PAGES 

26 

13.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

10.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

W.  SECURITY  CLASSIFICATION 

Or  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

Unlimited 

NSN  7540-01  280-5500 


Standard  Form  298  (rev.  2-89) 
Prescribed  by  ANSI  Std  Z39-18 
298-102 


Contents 


1.  INTRODUCTION  1 

2.  DESCRIPTION  OF  THE  LORAN  RADIOSONDE  SYSTEM  2 

3.  TEST  SERIES  4 

4.  DATA  ANALYSIS  5 

5.  results  ^  (pusnrfl®*010, 4 

6.  CONCLUSIONS  15 

APPENDIX  A:  SONDE/RADAR  STATISTICS  17 


Illustrations 

1.  Lor  an  Radiosonde  System  3 

2a.  urad  +  ugonde  vs  Time  (Before  Aligning)  7 

2b.  vrad  +  vs  Time  (Before  Aligning)  8 

3a.  urad  +  usonde  vs  Time  (After  Aligning)  9 

vrad  +  vsonde  vs  Time  (After  Aligning)  10 

4.  (urad  “  usonde)  +  (vrad  ”  vsonde) vs  Time  1 1 

5.  Histogram  of  (urad  -  u^n^)  12 

6.  Histogram  of  (vratj-vgo^g)  13 

7-  KadS  ~  urad3)  +  (vrad5  “  vrad3) vs  Time  14 

Tables 

1.  Vest  Summary  4 

2.  Results  15 

Al.  Sonde/Radar  Statistics  19 


iv 


Preface 


Frank  Schmidlin  from  NASA  Wallops  provided  guidance  In  setting  up  the  field  experiment, 
was  responsible  for  providing  the  range  tracking,  and  provided  specialized  processing  of  the 
radar  data,  which  facilitated  the  Intercomparison  with  the  radiosonde  data.  I  thank  him  for  his 
efforts. 


V 


Wind  Measurement  Accuracy 
of  a  Loran  Radiosonde 


1.  INTRODUCTION 


The  Atmospheric  Sciences  Division  of  Phillips  Laboratory  has  begun  a  program  to  see  If  wind 
measurements  made  by  a  wind  profiler  could  satisfy  day  of  launch  (DOL)  mlsslle-range  wind 
requirements.  The  wind  profiler  selected  for  this  was  one  of  the  403  MHz  production  units 
procured  by  NOAA  from  the  Unisys  Corporation. 

The  first  step  In  the  program  was  to  choose  and  evaluate  a  standard  against  which  the 
accuracy  of  the  wind  profiler  would  be  compared.  A  Loran  type  radiosonde  system  was  selected  as 
the  comparison  standard  for  several  reasons;  accuracy,  cost,  portability,  and  availability.  The 
results  of  the  test  series  In  which  the  accuracy  of  the  Loran  sonde  system  was  measured  Is  the 
subject  of  this  report. 

The  missile  range  requires  an  accuracy  of  1  m/s  with  a  resolution  of  about  200  m  so  it  is 
important  that  the  accuracy  of  the  testing  system  be  significantly  better  than  this.  The  accuracy 
of  the  Loran  system  Is  predicted  to  be  <  0.5  m/s  In  the  Sudbury  area  where  the  profiler  Is  to  be 
located  and  also  In  the  Wallops  Island  area. 

There  were  a  couple  of  Loran  radiosonde  systems  available  commercially  but  we  chose  to 
assemble  our  own  from  off  the  shelf  hardware.  This  allowed  us  to  customize  the  system  through 
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software  developmeni  to  satisfy  our  specific  requirement.  In  addition  it  allowed  us  to  build  two 
systems  with  an  outlay  of  about  a  third  the  cost  of  buying  two  complete  sonde  systems. 

The  Loran  systems  are  portable  in  that  they  do  not  rely  on  antenna  pointing  to  derive  winds 
as  does  radar.  Consequently,  a  Loran  system  does  not  need  a  large  antenna/pedestal  or  elaborate 
siting  procedures. 

Since  the  software  programs  used  to  track  the  Loran  signals  and  to  derive  winds  were  new.  it 
was  necessary  to  test  the  system  against  a  known  tracking  system  to  assess  its  accuracy.  This  was 
accomplished  at  Wallops  Island  using  Radar  5.  an  AN/FPQ-6  tracking  radar,  and  Radar  3.  an 
AN/FPS-16  tracking  radar,  in  June  1990. 


2.  DESCRIPTION  OF  THE  LORAN  RADIOSONDE  SYSTEM 


A  block  diagram  of  the  system  is  shown  in  Figure  1 .  The  component  parts  are  all  off  the  shelf 
or  were  available  in-house.  The  receiver  and  pre-amp  are  Communitronics  MR14  and  MR14A 
respectively,  which  were  available  in-house.  The  meteorological  processor  is  a  Vaisala  PTU 
Processor  PP1 1 .  The  Loran  tracker  is  the  tracking  cards  from  a  Viz  Manufacturing  Corporation 
W9000  system.  The  antenna  was  a  simple  corner  reflector  (=13  db  gain)  available  from  previous 
work.  The  system  uses  an  Arrowrotor  antenna  rotator  from  Radio  Shack.  The  PC/AT  was  a  286 
clone  from  the  government  desktop  contract. 

The  system  uses  a  Vaisala  RS-80-  15L  sonde  .  The  sonde  transmits  a  403  MHz  FM  modulated 
signal.  The  Loran  (100  kHz)  and  the  met  data  (8-10  kHz)  are  frequency  multiplexed.  The  met  data 
band  is  also  time  multiplexed  between  pressure,  temperature,  humidity,  and  some  calibration 
data.  The  receiver  demodulates  the  data  and  separates  the  met  data  and  the  Loran  data.  The  met 
data  is  provided  to  the  met  data  processor  which  converts  it  to  engineering  units  using 
calibration  data  from  a  paper  tape  provided  with  each  sonde.  The  met  data  processor  then 
outputs  this  data  as  an  RS  232  data  stream  to  the  PC/AT.  A  complete  frame  of  met  data  requires 
approximately  1 .2  seconds.  The  Loran  data  is  provided  by  the  receiver  to  the  Loran  tracker, 
which  converts  it  into  time  delays  (TD's)  and  provides  these  TD's  as  an  RS  232  data  stream  to  the 
PC/AT. 

The  PC/AT  ingests  Loran  TD's  every  2  seconds.  It  ingests  all  the  met  data  but  uses  only  the 
latest  met  data  available  every  2  s.  Both  the  2-s  met  and  the  2-s  Loran  data  are  written  to  files  on 
floppy  disks.  These  can  be  used  to  rerun  the  flights  with  different  parameters,  such  as  Loran 
stations,  smoothing  interval,  and  wind  calculation  interval.  The  system  also  uses  the  data  to 
generate  a  6-s  data  file  that  contains  time  of  flight  (s),  pressure  (mb),  alt  (km),  temp  (°C),  humidity 
(%RH),  X-distance  east  (km).  Y-distance  north  (km),  u-velocity  east  (m/s),  and  v-velocity  north 
(m/s).  These  6-s  data  are  written  to  the  floppy  disk  and  are  also  kept  in  active  memory  for 
plotting  during  flight. 

During  flight  the  operator  has  the  choice  of  two  types  of  displays.  One  has  two  plots  on  tire 
screen  together  with  a  data  window  that  continuously  displays  the  last  seven  6-s  data  frames  and 
a  menu  window.  The  two  plots  can  be  changed  by  the  operator  to  several  choices;  T.H  versus  alt. 
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u.v  versus  alt,  X  versus  Y.  P.T.H  versus  time,  X.Y  versus  time.  The  other  display  has  only  a  single 
plot  on  the  screen  which  can  be  any  of  the  plot  choices  given  above. 


3.  TEST  SERIES 


The  tests  were  performed  at  the  Wallops  Island  Test  Facility  on  5.  6.  and  8  June  1990.  Ten 
ascents  produced  useful  results.  Two  additional  ascents  were  not  used  due  to  Instrumentation 
difficulties.  Table  1  gives  details  of  the  ascents  used  In  the  data  analysis.  Ascents  1-8  were 
tracked  by  two  Loran  ground  stations.  Ascent  6  had  two  radiosondes  separated  horizontally  by 
an  8-ft  pole  and  offset  in  frequency  so  that  the  Loran  ground  stations  could  track  different 
sondes.  The  other  ascents  had  the  Loran  ground  stations  tracking  the  same  sonde.  Ascents  9  and 
10  had  single  Loran  ground  stations  due  to  a  malfunction  In  the  No.l  ground  station.  Ascents  8 
and  9  were  tracked  by  two  radars. 


Table  1.  Test  Summary 


Date 

Time 

Ascent  No. 

Sonde  No. 

Radar 

6/5/90 

1310 

1 

1 

i 

5 

1 

2 

5 

6/5/90 

1645 

2 

2 

1 

5 

2 

2 

5 

6/5/90 

1830 

3 

3 

1 

5 

3 

2 

5 

6/6/90 

1309 

4 

4 

1 

3 

4 

2 

3 

6/6/90 

1438 

5 

5 

1 

3 

5 

2 

3 

6/6/90 

1710 

6 

6a 

1 

5  double  sonde  on  bar 

6b 

2 

5 

6/6/90 

1858 

7 

7 

1 

5 

7 

2 

5 

6/8/90 

1601 

8 

8 

1 

3.5  double  radar  track 

8 

2 

3,5 

0/H/!X) 

1721 

9 

9 

2 

3.5  double  radar  lra<  k 

6/8/90 

1838 

10 

10 

2 

3 

4.  DATA  ANALYSIS 


The  wind  profiler  produces  a  wind  measurement  for  each  250  meter  Increment  of  height  from 
500  m  above  the  surface  to  16250  m.  These  winds  are  derived  from  radar  returns  from 
approximately  350  m  of  altitude  In  the  lower  beam  (500  m  to  9500  m)  and  from  approximately 
1050  m  of  altitude  In  the  upper  beam  (7500  to  16250  m).  I  decided  to  analyze  the  test  data  using  die 
difference  In  position  data  1  min  apart,  which  averages  the  sonde  motion  during  this  Interval. 
Since  the  radiosonde  rises  approximately  350  m/mln.  this  matches  the  altitude  Interval  of  the 
radar  returns  contributing  to  the  winds  In  the  lower  beam  of  the  profiler.  Modifying  this  liner 
scale  data  to  reflect  the  interval  In  the  upper  beam  is  simply  a  matter  of  vectorlally  averaging  3 
min  of  the  1-mln  data. 

The  radar  winds  were  determined  every  6  s  by  taking  the  difference  In  position  data  1  min 
apart.  The  position  data  used  was  the  result  of  averaging  5  s  of  1-s  data  (51  data  points). 


(1) 


(2) 

where  Xj.acj(l).  xra(j())  are  distances  east  of  the  launch  point  and  Yrad(l).  yratj(J)  are  distances 
north  of  the  launch  point.  Xj.ad(j)  an<*  Yrad^  316  °btalne<l  every  0. 1  s  and  X^^i).  Yrad<‘>  are 

obtained  every  6  s. 

Then 


25 


Xrad(i)  =  — ~  Xra<l(j) 


51 


j  =  -25 


,  25 

Yrad(i)  =  __  ^  yrad(J) 

51  |  =  25 


Uiad(l)  =  Xradfl  +  5)  -  Xm<j(i  ~  5) 


(3) 


vrad(l)  —  Yrad(i  +  5)  -  Yradfl  "  5)  . 


(4) 


The  radiosonde  winds  are  derived  In  a  similar  manner. 

|  =  5 

Xnonde(i)  =  ^  X»onde(J) 

)  =  -5 

1  =  5 

Yaondf(i)  =  ysonde(j)  • 

J  »-5 


(5) 

(6) 
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Since  the  radiosonde  raw  data  XgondeO)-  ysonde^  are  available  In  2-s  Intervals  this  represents  a 
20-s  smoothing  of  the  position  data  prior  to  differencing  for  the  winds.  The  wind  differencing  for 
the  radiosonde  Is  Identical  to  the  radar  Eqs.  (3)  and  (4). 

There  was  no  common  time  base  between  the  two  systems  so  It  was  necessary  to  shift  the  two 
wind  sequences  of  6-s  winds,  radar  and  sonde,  In  time.  This  was  accomplished  on  a  computer 
where  the  sequences  were  displayed  as  two  graphs  (time  vs  uratj(l)  and  usoncjeU)-  Ume  vs  vracj(l) 

and  vsoruje(i).  Figure  s  2a  and  2b  show  an  example  prior  to  shifting.  In  addition  the  statistics 

(Number  of  pts,  variance,  mean,  and  max  dlff)  were  presented  to  help  In  alignment  prior  to  the 
final  statistics  being  recorded.  Figures  3a  and  3b  show  the  same  example  after  shifting.  Figure  4 
shows  the  difference  between  the  radar  wind  values  and  the  radiosonde  wind  values.  The 
difference  in  the  u  components,  on  the  right,  have  been  given  a  +2  m/s  offset  while  the  v 
components  have  a  -2  m/s  offset,  to  avoid  confusion.  Figures  5  and  6  show  histograms  of  the  u 
and  v  differences  respectively.  Figure  7  shows  the  differences  between  the  two  radar  tracks  for 
comparison  with  the  differences  between  sonde  and  radar  Figure  4.  Figures  2  through  7  are  all  for 
ascent  8,  system  2. 


5.  RESULTS 


The  results  are  given  In  Table  2  as  rms  errors  in  m/s  for  each  system  and  each  ascent  for  both 
the  u  and  v  components  and  the  vector  sum  of  the  components.  In  addition  the  difference  of  the 
vector  rms's  from  system  1  and  system  2  are  given  where  appropriate.  There  are  five  flights  with 
a  single  sonde  tracked  by  both  Loran  ground  sets  and  radar  5.  The  average  vector  error  using 
system  1  Is  0.327  m/s  and  using  system  2,  0.333  m/s.  If  we  Include  all  ascents  iiji;ig  one  sonde  and 
two  Loran  ground  tracks,  regardless  of  radar,  we  have  seven  cases  and  an  rms  vector  error  of 
0.323  m/s  for  system  1  and  0.320  m/s  for  system  2. 

Using  all  the  ascents  (10)  and  the  system  2  data  the  rms  vector  error  Is  0.317  m/s.  The  average 
mean  difference  using  these  same  10  ascents  with  system  2  was  0.0953  m/s  In  the  u  axis  and 
0.0220  m/s  In  tile  v  axis.  Two  flights  were  tracked  by  both  radars.  The  average  vector  rms 
differences  between  the  two  radars  was  0.078  m/s  which  if  the  error  Is  ascribed  equally  would  be 
0.055  m/s  per  radar.  The  average  mean  difference  was  0.002  m/s  or  less  per  axis  for  these  radar 
data. 
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Figure  5. 
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Figure  6.  Histogram  of  (vracj  - 


Table  2.  Results 


System  1 

System  2 

Ascent 

Number 

Radar 

Number 

u 

rms 

V 

rms 

Vector 

rms 

Diff 

Vector  rms 

u 

rms 

V 

rms 

Vector 

rms 

1 

5 

0.220 

0.186 

0.288 

0.061 

0.248 

0.245 

0.349 

2 

5 

0.286 

0.221 

0.361 

0.003 

0.300 

0.206 

0.364 

3 

5 

0.278 

0.256 

0.378 

-.022 

0.278 

0.222 

0.356 

4 

3 

0.249 

0.238 

0.344 

-.035 

0.248 

0.184 

0.309 

5 

3 

0.232 

0.159 

0.281 

-.013 

0.224 

0.147 

0.268 

6* 

5 

0.204 

0.159 

0.259 

0.076 

0.268 

0.201 

0.335 

7 

5 

0.201 

0.194 

0.279 

0.008 

0.226 

0.177 

0  287 

8 

5 

0.279 

0.172 

0.328 

-.020 

0.252 

0.177 

0.308 

8 

3 

0.282 

0.174 

0.331 

-.020 

0.256 

0.177 

0.311 

9 

5 

0.302 

0  180 

0.352 

9 

3 

0.298 

0.180 

0.348 

10 

3 

0.273 

0.209 

0.344 

*  Separate  sondes  for  system  1  and  system  2. 


The  appendix  contains  a  more  complete  listing  of  statistics  for  each  flight  on  which  this  table 
is  based. 


6.  CONCLUSIONS 


These  tests  show  an  overall  rms  vector  difference  of  0.32  m/s  between  winds  obtained  using  a 
Loran  system  and  winds  derived  using  tracking  radars  at  GSFC/Wallops  Flight  Facility.  The 
difference  between  winds  derived  using  the  two  radars  tracking  the  same  sonde  results  in  an  rms 
vector  difference  of  0.055  m/s  per  radar.  This  radar  uncertainty  accounts  for  only  0.005  m/s  of 
the  0.32  m/s  difference  between  the  sonde  and  radar.  Both  the  radar  and  sonde  winds  represent 
position  differences  1  minute  apart.  This  represents  approximately  350  meter  vertical 
resolution.  The  radar  position  data  used  in  these  calculations  were  generated  by  averaging  5  s  of 
data  while  the  sonde  position  data  represents  20-s  averages. 

There  were  two  ground  stations  tracking  a  single  sonde  and  there  appears  to  be  no  significant 
difference  due  to  the  ground  station  used  (0.328  vs  0.323  m/s).  There  was  only  one  ascent  where 
different  sondes  were  tracked  by  each  ground  station.  The  difference  (0.076  m/s)  in  the  rms  vector 


differences  was  much  larger  than  the  differences  when  tracking  the  same  sonde,  but  the 
magnitude  is  unimportant  in  terms  of  radiosonde  accuracies. 

The  accuracy  of  the  radiosonde  (0.32  m/s)  is  sufficient  to  use  for  evaluating  the  profiler  since 
its  variance  (0.102  m2/s2)  would  be  only  l/10th  of  the  variance  (1.0  m2/s2)  budget  needed  to 
demonstrate  whether  the  profiler  can  meet  a  1 .0  m/s  requirement. 
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Appendix  A 


Sonde/Radar  Statistics 
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